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Some numerical studies of the evolution 
of generalized parton distributions 
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Abstract: We study the evolution behavior of generahzed parton distributions at smaU longitudinal mo- 
nicntuni fraction. Particular attention is paid to the ratio of a generalized parton distribution and its forward 
limit, to the mixing between quarks and gluons, and to the dependence on the squared momentum transfer t. 



1 Introduction 

A characteristic property of generalized parton distri- 
butions (GPDs) is their renormalization scale depen- 
dence, described by evolution equations whose deriva- 
tion led to the very discovery of these functions more 
than a decade ago [1]. On a practical level, the scale 
dependence of GPDs is of direct importance for the 
quantitative description of exclusive scattering pro- 
cesses. Moreover, understanding general features of 
the evolution behavior should be helpful for developing 
realistic models and parameterizations of GPDs. The 
question how a given input distribution changes when 
evolved to higher scales has been addressed in several 
studies, both numerically and analytically [2-8]. Fur- 
ther progress has been achieved recently [9,10] by con- 
structing explicit solutions of the evolution equations 
with methods that generalize the familiar Mellin mo- 
ment inversion for parton density functions (PDFs). 

The aim of the present contribution is to study a 
number of aspects in the evolution of GPDs at a nu- 
merical level. We will largely concentrate on the value 
of the GPDs at a; = ^, which at leading order in as de- 
termines the imaginary part of scattering amplitudes, 
and via dispersion relations also gives their real part 
up to a ^ independent constant [11]. Furthermore we 
will focus on the region of small i^, where the behavior 
of distributions can be conveniently approximated by 
a power-law behavior. We will pay special attention 
to the mixing between the gluon GPD H^{x,£,,t) and 

"/ 

H'^ix, tt)=Y. t^'^^' ^' ^) - ^''(-^' ^' *)] ' (1) 
q 

whose forward limit 

"/ 
H'^ix, 0, 0) = Six) = J2 [«(^) + «»] (2) 



is the familiar singlet combination of quark and an- 
tiquark PDFs. For comparison we will also consider 
H''-'^{x,^,t) = H''{x,^,t) - H'^{x,£,,t) as a represen- 
tative of the non-singlet sector. 

After specifying in Sect. [2] the GPD model used as 
initial condition for the evolution, we devote most of 
Sect. [3] to a quantitative study of the old question how 
the ratio of GPDs and PDFs behaves when evolved 
to higher scales. We shall in addition take a look at 
the behavior of the GPDs around x — ^. In Sects. |4] 
and [5] we turn to the dependence of GPDs on the 
squared momentum transfer t. Both theoretical con- 
siderations [12] and lattice QCD calculations [13] in- 
dicate that this dependence is correlated with the one 
on the longitudinal variables x and ^. Since evolution 
affects the x dependence at given ^ and t, it also af- 
fects the t dependence at given x and ^ in a nontrivial 
fashion, which we will quantify in two model scenarios. 

For our calculations we have used the numerical code 
of [14], which provides a numerically fast and sta- 
ble implementation of GPD evolution at leading or- 
der (LO) in ag. The effects of next-to- leading (NLO) 
and next-to-next-to-leading order (NNLO) terms in 
the evolution kernels have been studied [8,15] and are 
known to be important, especially at small ^ in the 
gluon and singlet sector. This should be kept in mind 
as a caveat when interpreting our results, but we think 
that a study at LO is still of some relevance. On one 
hand, the arguments in [5,7] about the pattern of evo- 
lution to higher scales are based on the LO kernels, so 
that this order is adequate to test the numerical valid- 
ity of these arguments. On the other hand, evolution 
effects on the t dependence are barely known at all, 
and LO results should at least provide a valid starting 
point for further investigation. 

In the evolution kernels and the running coupling 
we take Uf — A for uic < /i < mi, and Uf = 5 for 



M ^ "^fci with the charm and bottom quark masses 
mc = 1.3 GeV and mt = 4.5 GeV used in the CTEQ6 
parton analysis [16], which we use for calculating the 
model GPDs at the starting scale of evolution. We fur- 
thermore follow the CTEQ6 analysis in taking the two- 
loop running coupling with A*^^^ = 326 MeV and A^^^^ — 
226 MeV, which corresponds to a'^^rric) = 0.40 and 
a'-g^^AIz) — 0.118. We shall not consider scales be- 
low fi — 1.3 GeV, which we regard as a compromise 
between starting evolution at a "low scale" and stay- 
ing in a region where as is not so large that the LO 
approximation becomes more and more questionable. 

2 Initial conditions 

At the starting scale of evolution, we use the Musatov- 
Radyushkin ansatz, which is based on double distribu- 
tions [7] . With the conventional definitions of H'^ and 
H^ , given e.g. in [17], we can write this ansatz as 



/•I /-i-l/^l 

m{x,S„t) ^ dp da 

J-l "'-1 + 1/31 



S{x-P-^a) 



Ha{x,i,t) 



X /i6(/3,a)iJ«(/3,0,i) 



1 /-i-l/SI 
dp / da5{x-P- ^a) 

-1 "'-1 + 1/31 



xhbiP,a)Ha{f3,0,t) 



with 



hb{P,a) 



r(26 + 2) [{i-\p\r 



„2lb 



22fc+ir2(5-f 1) (l-|/3|)2b+l 



(3) 



(4) 



In this work we will use different values of the profile 
parameter b, which for simplicity will always be taken 
equal for all quark and gluon distributions. The ansatz 
^ has been extensively used in the literature so far. 
One should keep in mind that it does not exhaust the 
possibihties of modehng, and other approaches [18-21] 
are being pursued in the literature. As we will see, 
this model does however provide enough flexibility to 
address a number of important questions. 

The model also permits useful analytic approxima- 
tions at small ^. At a; = ^ the integrals in Q are 
restricted to /3 < 2^, so that for ^ <C 1 one can neglect 
the P dependence in /ib(/3, a). Approximating 1 + ^ by 
1 in the integration limits, one then has [7] 



H'{^,^,t) 



jj dphJo,l-^)H\P,0,t) 
r(26 + 2) '■1 



T^ib + 1) Jo 



dz{l-z)''z''H\2^z,0,t) (5) 



with i = q,g. We will use this approximation shortly. 

3 Evolution at fixed t 

In this section we study the evolution of GPDs at a 
fixed value of t. We take t — and do not display this 
variable for brevity. To quantify the difference between 
generalized and usual parton distributions we use the 
conventional skewness ratios 



R^^,^i) 



i?'(e,Ai) 






(6) 



where we have explicitly displayed the dependence on 
the scale /x in the distributions. 

As is well known, the PDFs obtained from fits to 
data follow an approximate power-law behavior at 
small X, 



xg{x) 



„-A 



xq{x) 



-A 



(7) 



at given /i, where a and A depend of course on the 
parton species. With the ansatz ([3]) for GPDs this 
leads to a power-law behavior 



ff^(c,o^r 



^H^i^O-C 



(8) 



of the GPDs at small ^ according to ([5]), with the same 
powers A as for the corresponding PDFs. The skewness 
ratios at small ^ are readily obtained as 



RlW 



^g^,,_ r(26 + 2) r(& + i-A) 

''^"' T{2h+l-\) T{h+1) 



^^"' T{2b + 2-X) r(5+l) 
r(26 + 2) r(6-A) 



R1W 



(9) 



at the scale /x where the ansatz ([3]) is made. Numeri- 
cally, we find that the approximate power-laws ([7]) and 
([5|) remain valid under evolution to higher scales, with 
powers A that depend on fi but remain the same for 
the forward distributions and the GPDs. 

Based on the considerations using the Shuvaev 
transformation, it has been argued in [5] that at small 
f and high enough scale, the skewness ratio should be 
given by [4, 22] 






(A) = 



22A+3 r(A + f ) 



RkW 



V^ r(A + 4) ' 

22A+3 r(A + f) 



/¥ r(A + 3) 



(10) 



for gluons and quarks, respectively. Here A is the 
power in ([7]) at the scale where i?^(^,/i) or i?*(^,/i) 



is evaluated. More precisely, the ratios in (|T0|) are ob- 
tained if ([7]) holds and if all Gegenbauer moments of 
the GPD in question are independent of £,. Musatov 
and Radyushkin [7] have shown that at small x and ^ 
this condition is tantamount to the GPD being given 
by ([2]) with 6 = A + 1, for both gluon and quark dis- 
tributions. Indeed, one can easily check that 



^Sh('^) - ^A+l('^)- 



(11) 



for i — g,q. Using a different line of arguments, the 
authors of [10] also expect that pO)) should become 
valid after LO evolution to high scales, provided that 
one takes a particular joint limit of large fi and 1/^. 
The relations (fTO|) are often used to calculate high- 
energy scattering amplitudes, so that it is important 
to test under which conditions they may be assumed 
to hold. 

We have taken the double distribution model ([3]) 
with the CTEQ6L distributions [16] at fiQ = 1.3 GeV 
as input. After LO evolution to a scale /i, we have fit- 
ted effective power laws for g{x) and i/^(^,^), and we 
have evaluated the skewness ratio R^{^,fi) from ^. 
In analogy we have determined power laws and ratios 
R^ and i?""'' for the combinations H^ and H"^'^ in- 
troduced in Sect.[T] 

Let us first discuss the power-law behavior ([7]) of 
the PDFs, which is not exact and only valid in a 
certain range of x. We fitted power-laws to the 
CTEQ6L parameterization for g{x), S{x) and u{x) — 
d{x) in the three intervals [10-^10-*], [lO""^, 10"^] 
and [10~^,10~^]. The resulting powers for the gluon 
and quark singlet distributions are shown in Fig. [TJ 
We see a clear x dependence of the effective power A, 
especially at larger x. In Fig. [3 we show the powers 
obtained in the interval 10"'* < x < 10"^ for a larger 
range of /i. We note that under evolution the powers 
for the gluon and the quark singlet become similar but 
remain different up to very high /i. This effect has al- 
ready been pointed out in [23] . For the non-singlet dis- 
tribution u ^ d the effective power A is between —0.41 
and —0.42 in all three x intervals. It changes by less 
than 1% under evolution in the fi range corresponding 
to Figs. [Hand [21 

According to (O there is no simple relation between 
the ranges of x and ^ in which the same power-law 
behavior should approximately hold for a PDF and 
the corresponding GPD. For simplicity we have fitted 
iJ9(C,C), -ff^(C,0 and H^'-'^iC^O 1° power laws (0) in 
the same ^ intervals that we took for the PDFs. An 
example of such a fit is shown in Fig. [3l where we see 
that H^{^, ^) indeed follows an approximate power-law 
over about one order of magnitude in ^ but not over 
the full range of the plot. We find that corresponding 
powers A for PDFs and GPDs differ by at most 3% in 



the respective x and /i ranges of Figs. [T] and [2l An 
exception is the quark singlet distribution in the inter- 
val 10"^ < a; < 10~^, where the power for the GPD 
is higher than that for the PDF by 5% to 10%. This 
is not surprising, given that already in Fig. [T] we see 
a more rapid change of the effective power at higher 
X. For definiteness we will evaluate Rb{X) and i?sh(A) 
with the powers fitted to the PDFs. We have checked 
that our conclusions do not change when taking the 
powers for the GPDs instead. 

We note that in a specific joint limit of large fi 
and 1/^, the solutions of the LO evolution equations 
for PDFs exhibit so-called double logarithmic scal- 
ing [24]. In this case one obtains d\n{xg)/di ~ i~-^/^ 
and dln{xS)/de - dln{xg)/d£ « -(2€)-i, where £ = 
lii^xo/x) with some constant xq. The effective pow- 
ers in (O are then larger for the gluon than for the 
quark singlet distribution and depend logarithmically 
on X. Double logarithmic scaling for GPDs in the re- 
gion X > ^ has been discussed in [10]. 

We have evaluated the skewness ratios R{^, /i) from 
the evolved GPDs and PDFs for ^ = 3.2 x 10"^ 
3.2x10"'' and 3.2xl0~3. This is compared with Rb (A) 
and -Rsh(A) calculated with A from our fits of the PDFs 
at the corresponding scale fi and in the corresponding 
X interval [10-^10-''], [10-^,10-3] and [lO^^^lO-^]. 
The result for gluons is shown in Figs. |4| and [S] where 
in the initial condition we have taken b — 2. At the 
starting scale ^o = 1.3 GeV the curves for i?^(^,/i) 
and R2{^) coincide as they should, whereas for increas- 
ing fj, they become different. This means that one ob- 
tains different results for if ^(^, £,; /i) when making the 
ansatz (|3]) at scale fi or when making it at scale /iq 
and then evolving the GPD. The difference is however 
fairly small. 

The curves for i?3(^,/i) and those for i?|i,(A) = 
R\+iW in Figs. [Hand [Dare rather close to each other. 
At the starting scale they hardly differ at all, which re- 
flects a particularity of the model ansatz Q for gluons. 
This is because the ratio R^{X) in ^ has a very weak 
b dependence for small A: varying b from 1 to cx) one 
obtains for instance _R^(0.1) between 1.072 and 1.088. 
With increasing A the b dependence grows only slowly, 
with i?^(0.3) between 1.231 and 1.307 for b between 1 
and oo. This is also seen in the left panel of Fig. [5l 
where i?^(A) is given for several b values. The solid 
curve in this figure is for 6 = 2 in the initial condition, 
but the corresponding results for 6 = 1 or & = 8 differ 
by less than 0.5% in the /i range of the figure. 

Obviously it is hard to see whether R^{£,, /i) tends to 
_Rgjj(A) under evolution if the two functions are already 
close at the starting scale. To investigate this further 
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Figure 1: Effective powers A obtained from fits xg{x) ^ x~^ (left) and xS{x) 
of X at given scale /n. 
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Figure 2: Effective powers as in Fig. (TJ obtained from 
fits in the interval 10~^ <x < 10"'^. 
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Figure 3: Values of Hs{£^,£^) (points) and power- 
law fits (lines) in successive intervals [10~^,10~^], 
[10-"*, 10-3], [10-3, 10-2] of ^. The lower curve is for 
the starting scale /Iq — 1.69 GeV^ and the upper one 
for ^2 = 50GeV^ 



we take a variant of ^, namely 

HS{x,^,t)^x dp daS{x~f3~^a) 

J-l "'-1 + 1/31 



X hbiP,a)f3~^H3{f3,0,t) 



(12) 



with hb{(3,a) as in @. This corresponds to a double 
distribution representation for x^^H^(x,^,t) instead 
of H^{x,£,,t), and one readily verifies that it gives 
Mellin moments of H^{x,£,,t) with a polynomial de- 
pendence on ^ as required by Lorentz invariance. An 
analogous representation was first discussed for the 
quark GPD of the pion [25] and was recently found 
to be relevant for polarized gluon GPDs [26]. In the 
case of H^ the ansatz (IT^ has the peculiar property of 
giving a zero at a; = that is not required by symme- 
try and quickly disappears under evolution. One may 
therefore not take this model too seriously, but it serves 
the purpose of giving a skewness ratio sufficiently dif- 
ferent from the one obtained with the more conven- 
tional ansatz ^. This is shown in the right panel of 
Fig. [51 where the dot-dashed curve corresponds to ini- 
tial conditions ^ for H^ and © for H'^ , with 6 = 2 
in both cases. We see that the ratio R^ in the two 
models indeed tends to a common value after evolu- 
tion. This value it not exactly equal to -Rgij(A) but 
differs from it by less than 2%. Such a small difference 
should not be regarded as significant: the form (fTU)l of 



i?gjj(A) is obtained in [4,22] from an integral of g{x) 
over X from ^/2 to 1, assuming the power behavior ([7]) 
in the entire interval. This is clearly an approximation. 
We now turn to the skewness ratio for the quark sin- 
glet distribution, which is shown in Fig. [5] In contrast 
to the gluon case, different values of b in the ansatz 
^ lead to significantly different skewness ratios at 
the starting scale. Evolution to higher /i brings the 
curves oi R^{^,fi) for different initial conditions closer 
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Figure 4: The skewness ratio for gluons. Here and in Figs.[5l[6]and[8l solid lines give the ratio i?(^, ^) of GPDs 
and PDFs evolved to the scale /x, and dashed or dotted lines give Rb{X) calculated with the effective power A 
fitted at that scale. 
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Figure 5: The skewness ratio for gluons at ^ = 3.2 x 10^''. The dot-dashed curve in the right plot is for the 
alternative model (fT2|) described in the text. 
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Figure 6: The skewness ratio for the quark singlet distribution. From top to bottom the solid curves correspond 
to 6 = 1, 2, 8 in the initial condition ^. At the starting scale \iq = 1.3 GeV they coincide with the corresponding 
dashed curves for i?f (A). 
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Figure 7: The quark singlet GPD for ^ = 3.2 x 10~* at the starting scale /Xg = 1.69 GeV^ and after evolution 
to /x^ = lOGeV . For each scale the upper curve gives H^{x,S,) and the lower curve gives {d/dx)lnH^{x,£^). 
The left plot is for 6 = 1 and the right plot for 6 = 2 in the initial condition. 



to each other. As in the case of gluons, they do not ex- 
actly approach the curve we calculate for i?fi,(A), but 
again this should not be regarded as significant since 
the power-law ([7]) with a fixed value of A is only an 
approximation for a certain x range. 

In the left panel of Fig. [S] we also see a clear difference 
between the evolved ratios R^{£,, IJ.) and the curves for 
R^{X) with A taken at the corresponding scale /i. In 
general there is hence a notable dependence of H^ on 
the scale where the ansatz (jS]) is made, especially for 
larger values of b. We find the dependence less pro- 
nounced for 6 = 1, in agreement with what was found 
in [7] for C = 5.26 x IQ-^. 

Based on the inversion of Gegenbauer moments, it 
was argued in [10] that quark distributions H'^{x,^) 
should develop a singular derivative {d/dx)H'^{x, ^) at 
a; = ^ after evolution. To investigate this, we have 
numerically calculated {d/dx)H^{x, £,) from the differ- 
ence quotient for successive points in x, which around 
X = ^ were spaced in intervals of 9 x 10^^. In Fig. [7] we 
plot H^{x,S,) together with its logarithmic derivative 
{d/dx)ln.H^{x,S,). Taking the derivative of ([3]) and 
making the same approximations which lead to ([5]), 
one finds that {d/dx)H^{x,£^) is singular at a; = ^ for 
b < 1 + A. Indeed, we see in the figure that at the start- 
ing scale the derivative has a singularity for 6=1 but 
remains finite for b — 2. Under evolution a singularity 
develops for 6 = 2, whereas for 6 = 1 the logarithmic 
derivative hardly changes. Notice that in the curves 
for H^{x, ^) one can barely recognize that the tangent 
at X = ^ should be vertical: this illustrates the limi- 
tations of rendering a weakly singular derivative in a 
plot. In contrast to the quark case, the ansatz ([3]) with 
6 = 1 or 6 = 2 gives a finite value of {d/dx)H^{x,£,) 
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Figure 8: The skewness ratio for the non-singlet distri- 
bution H"^'^. The meaning of the curves is as in the 
right panel of Fig. [51 



at x = ^. We have checked numerically that for both 
initial conditions the derivative remains finite under 
evolution, in agreement with what one expects from 
the analytical representation in [10]. 

To conclude this section, we briefly investigate the 
quark non-singlet distribution H"~'^. In Fig.[5]we show 
the skewness ratio for different values of 6 in the initial 
condition. The corresponding curves for R^~ (A) are 
not shown: they coincide with those for i?"^'*(^, n) at 
the starting scale and then remain essentially flat since 
the effective power A hardly changes with /i in this case. 
Under evolution to high scales the curves for different 6 
approach each other and the one for R^^ ( A) , although 
much more slowly than for R^ or R . 



6 



4 Ansatz for the t dependence 

To investigate the change of the t dependence with 
evolution, we wiU use the model ^ with 6 = 2 for 
quarks and gluons. We thus need an ansatz for the 
GPDs at zero skewness £, but finite t, which is described 
in this section. In all cases we assume an exponential t 
dependence that is correlated with x. For the valence- 
type combination of GPDs we take the form proposed 
inRef. [27]: 

H^{x, 0, t) = iJ«(x, 0, t) + iJ«(-x, 0, t) 

= qv{x) exp[tfq{x)] (13) 

with x > 0, qv{x) — q{x) — q{x) and 
f,{x)^aUl-xf\n- 

X 
+ Bq{l~xf+AqX{l~x)^. (14) 

The values < = 0.9 GeV"^ B„ = Bd = 0.59 GeV"^ 
An = 1.22 GeV"^ and Ad = 2.59 GeV~^ together with 
the CTEQ6M parameterization for qv{x) at /i = 2 GeV 
lead to a good description of the data for the Dirac 
form factors Fi (t) of proton and neutron, which are ob- 
tained by combining / dx H^{x, 0, t) for u and d quarks 
with the appropriate charge factors. 

For small x we can approximate p^ as fq{x) w 
a'^ ln(l/x) -|- Bq and thus have 



with 



where 



fqiO=a'Jn- + Bq, (19) 



Bq = Bq + a', (j2 ^^^ - l^^V (20) 



Turning to the gluon distribution, we take 

HS{x,0,t) = xgix) exp[tfgix)] (21) 

with the function 

fg{x)^a'g{l-x)Hn^ + Bg{l-x)\ (22) 

which has one parameter less than its counterpart (|14p . 
Since most phenomenological information about glu- 
ons presently comes from small-x data, it would be 
difficult to constrain a third parameter. The analog of 
(fT6l) reads 



(23) 



iine,e,i)~2^g(20exp[t/,(20] 
60 



(3 - A - to' )(4 - A - to' )(5 - A - to' ) 



for 6 = 2 and was already used in [28]. With the 
approximation in (|17p we find 



H^{x,0,t)^q4x)x~'<e'''- 

^ -{l+X+ta'J tBg 



(15) 



where in the second step we have assumed a small- a: 
behavior of the valence quark distributions as in ([7]). 
Since the x dependence of (fTSJ) is a power-law, the 
integral in ([5]) can be performed as in Sect. [3l and 
we can use ([9]) for the skewness ratio at small ^ after 
replacing A with A -I- ta'^ . For 6 = 2 this gives 



where 



and 



Ha{^,tt)^H\^,^,0)exp[tfgiO] 



m)=a'gln-+Bg 



Bg = Bg + a 



9 Z^ „ _ 



^n— 3 



n — A 



In 2 



(24) 



(25) 



(26) 



H^it^,t)^q,i20exp[tfqi20] 
60 



in analogy to the quark case. For our numerical study 
we take the parameters a' — 0.164 GeV^^ and Bg — 
1.2GeV~^ in order to match recent HI data on J/$ 
photoproduction, whose t dependence is well fitted by 



(2 - A - to;)(3 - A - to;)(4 - A - to^) 
For small t we can write 



(16) 



da 

—— ex exp 

dt 



6o + 4a;in^ 



(27) 



with values 6o = 4.63 GeV^^ and a'g = 0.164 GeV" 



1 



1 



n — \ — ta' n — X 
1 



exp 
exp 


-ln(l- 

K-x\ 


ta' \' 


n-XJ 



(17) 



n — A 
and thus approximate (I16p by 

mit^,t) « ff«(e,e,0) exp[i/,(C)] (18) 



for Wo = 90 GeV [29]. To connect 1^ with 
we have used the approximate relation da /dt oc 
|i7^(^,^, t)!'^, which is obtained at tree level when 
one keeps only the imaginary part of the scatter- 
ing amplitude. The skewness variable is given by 
2^ = [Mj/^i, /W^pY in terms of the ^p cm. energy. For 
simplicity, we have omitted the terms with l/{n — A) 
in ([26|l when fixing Bg. For typical values of A they 
are quite small. 
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For antiquarks we set 



H'i{-x,0,t) = -qix)exp[tf^ix)] 



(28) 



with X > 0. Little is known to date about the t de- 
pendence in the sea quark sector. Constraints can be 
provided by deeply virtual Compton scattering [30,31], 
which at small x is sensitive to both sea quark and 
gluon distributions. A comprehensive analysis of this 
data, as has recently been performed in [20], is be- 
yond the scope of this work. We will instead explore 
the pattern of evolution for two extreme choices. In 
model 1 we set the t slope fq equal to the one for va- 
lence quarks: 



Ju Ju 1 



Id Id 5 



Js Id 5 



(29) 



where the choice /j — fd has no strong motivation, but 
does not strongly influence the results we will obtain. 
In model 2 we set instead 



Jq — Jg 



(30) 



for all quark flavors. The initial conditions for evo- 
lution of the singlet and non-singlet combinations are 
then obtained from 

iJ«(x,0,i)= Y. [H^{x,0,t)-2H\~x,0,t)], 

q—u.d,s 

il"-^(x, 0, t) == H^{x, 0, t) - ff"(-x, 0, t) 

-H^{x,0,t) + H'^{-x,0,t). (31) 

For the evolution study in the next section, we make 
the ansatz ^ with the CTEQ6M parton distributions 
at /io = 2 GeV, so that we can use the fit of [27] for 
the t dependence of H^{x,0,t) as specified at the be- 
ginning of this section. In (f3T|) we have neglected the 
tiny charm quark distribution at fio- To explore the 
region of lower scales, we will also consider backward 
evolution. 

5 Evolution of the t dependence 

In accordance with the analytical considerations in the 
previous section, we find that at the initial scale the 
t dependence of H3{£^, ^, t) is well described by an ex- 
ponential form at small t and ^. Evolving to higher 
scales we still find an approximately exponential be- 
havior in both model 1 and 2, as shown in Fig. [9] for 
^ = 3.2 X 10^^. A slight departure from an exact 
exponential in the full region < — i < 1 GeV is 
however visible at fp = 50GeV^. Evolving to lower 
scales, we still find an approximate exponential t de- 
pendence at /z^ = 3 GeV , but for yet lower scales the 
situation changes. At fi^ — 2 GeV" the distribution 



H^{£,,£,,t) turns negative for —t around 0.5 GeV in 
model 1 and around 0.3 GeV in model 2, whereas at 
fi^ = 1.69 GeV^ we have Hs{^,^,t) < already for 
t = 0. This is due to the behavior of the CTEQ6M 
gluon density at low scales. Since the gluon distribu- 
tion in this region varies considerably between different 
global parton fits, we shall not elaborate on this issue 
further here. 

The singlet distribution H^ {(_,£,, t) is again well ap- 
proximated by an exponential in t at the starting scale, 
and it stays exponential to high accuracy in model 2 
up to fi"^ = 50 GeV^ and even down to fj,"^ — 2 GeV^ . 
As shown in Fig. [9l this is however not the case in 
model 1. Here we find a clear departure from an expo- 
nential behavior even when evolving from the starting 
scale to /i^ = 6 GeV^ , whereas under backward evolu- 
tion _ff (^,^,t) rapidly turns negative for some value 
of t. We notice that in model 1 the x dependence of 
H^{x,0,t) at the starting scale rapidly changes with 
t due to the large value of a'^. This induces a corre- 
sponding change in the x dependence of H {x,£^,t), 
which enters in the evolution equations. 

To quantify the change of the t dependence under 
evolution, we fit the GPDs at given ^ and /i to 

Wi^,^,t■,^i)^W{^,tO■,^l) exp[i/,(e;M)] (32) 

for —t between and 0.5 GeV^, where i — g, S. Given 
the behavior of the distributions under backward evo- 
lution, we restrict these fits to /i^ > 4GeV^. Whereas 
for H^ in model 2 and for H^ in both models the 
form (j32p gives an excellent description in the kine- 
matical region of the fit, the corresponding fit for H 
in model 1 can only be approximate, as is seen in Fig.[9l 
This must be kept in mind when interpreting the sub- 
sequent results, but despite this caveat the correspond- 
ing t slope fs gives a fair account of how H^ {£_,£_, t) 
changes with /i. The results of the fit are shown in 
Fig-Hnifor the starting scale and for ^^ — 50 GeV^. We 
see that over a wide region of small ^ the dependence 
of fi{(.', A*) on ^ remains logarithmic after evolution to 
higher scales. For given fi we can hence perform a fit 

M^;^i)^a'M\n^+B,{f,). (33) 

The results of such a fit in the range 3.2 x 10~^ < 
^ < 3.2 X 10~^ are shown in Fig. [Til where we plot 
/i(^;/^) at the midpoint ^ = lO^"' of the fit range, 
as well as the effective shrinkage parameter a^(/i). In 
model 2, fi{^; fJ.) and a^(/i) are equal for the gluon and 
the quark singlet to a good precision at the starting 
scale by construction. They change rather mildly un- 
der evolution to higher scales, but a visible difference 
between gluon and singlet appears, especially for a^. 
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Figure 9: H^{£_,£^,t) and H^{£_,S^,t) at ^ = 3.2 x 10^** for different scales fi^ as indicated. Here and in the 
following figures, solid curves are for model 1 and dashed curves for model 2. At the starting scale ^^ = 4 GeV^ 
the curves for H^{^,£^,t) coincide in both models. 
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Figure 10: The t slope /(^;/i) fitted according to ([32|) for the gluon (left) and the quark singlet GPD (right). 
At /i^ ~ 4GeV^ the curves for fg{S,',l^) coincide in models 1 and 2. 
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Figure 11: The t slope /i(f ; /i) at ^ = 10 ^ and the shrinkage parameter a'^ in (|33p as a function of /i^ 
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Figure 12: iJ"-''(^,^, t) at ^ == 3.2 x 10"'* for different 
scales fi^. 



In model 1, we see that the slope and the shrinkage pa- 
rameter for the singlet evolve quite strongly and tend 
to approach the corresponding values in the gluon dis- 
tribution, which increasingly dominates evolution with 
increasing /x. The respective values for the gluon and 
the quark singlet are however clearly different even at 
^2 = 50GeV2. 

We have also investigated the evolution behavior of 
the non-singlet quantity H^~'^{^, ^, t). Only the flavor 
difference enters for the sea quark distributions ([28]) in 
this case, and we restrict our investigation to model 1. 
The t dependence changes in a similar way as for the 
quark singlet in model 1 , as becomes evident from com- 
parison of Fig. [S] and Fig. [T^ In particular, evolution 
to higher scales modifies the exponential behavior of 
the initial condition. A fit of the t slope for —t be- 
tween and 0.5 GeV must hence be taken with the 
same caveat as above. The t slope fu-d{^', fJ-) fitted as 
in p2p shows an approximately logarithmic ^ behavior 
in the fuh range 2GeV^ < fJ-'^ < 50GeV^, so that we 
can again perform a fit to the form p3|) in the region 
3.2 X 10"^ < ^ < 3.2 X 10"''. Between ^^ ^ 2 GeV^ 
and 50 GeV , the resulting shrinkage parameter a'^_^ 
increases only by about 2%, and the slope /u-d(C;M) 
at ^ = 10"** decreases by about 10%. Compared with 
the quark singlet sector in model 1, evolution effects in 
the non-singlet sector are hence considerably weaker. 



6 Conclusions 

We have studied several aspects of the evolution be- 
havior of GPDs at small £,. To do this, we assumed 
a particular form of the GPDs at a moderately low 
scale and numerically evolved this model ansatz to 
higher fi. At i = we have taken initial conditions 
for which i?(^, ^, 0) and the corresponding parton den- 



sity H(^, 0, 0) approximately obey power-laws with the 
same power. Under evolution to higher scales this 
power changes but remains the same for a GPD and 
the associated PDF. As a consequence, the skewness 
ratio R{^, ^i) is only weakly ^ dependent, to the extent 
that the effective power changes with f . The values 
of R{S,, /i) for different initial conditions approach each 
other with increasing /i, and at high scales they are 
well approximated by the Shuvaev formula (jlOp . This 
convergence is however not very fast: with rather dif- 
ferent values of i?(^, //) at /i^ = to^ it only becomes 
visible at /i^ of a few 10 GeV^ for the gluon and quark 
singlet distributions, and at yet larger values in the 
non-singlet sector. We have not attempted to study 
how the situation would change for initial conditions 
at much lower scale, considering that in this case the 
leading-order approximation of the evolution equations 
would no longer be suitable for drawing quantitative 
conclusions. We confirm the finding of [10] that evo- 
lution to higher scales generates a singular derivative 
{d/dx)H{x,£,,t) at a; = ^ for quarks, but not for glu- 



To study the change of t dependence under evo- 
lution, we have chosen initial conditions at /io — 
2 GeV such that H{C,tt) -^ exp[</(^)] and /(O = 
a' ln(l/^) + i? at small ^ and t. For distributions 
with a small shrinkage parameter a' , we find that to 
a good approximation the t dependence remains expo- 
nential under evolution to higher (and to some extent 
also to lower) scales. In contrast, a deviation from 
an exponential t behavior becomes visible after evolu- 
tion rather quickly for distributions with large a' at 
the starting scale, so that a fit to an exponential form 
is only approximate in these cases. The fitted slopes 
f{£_) of the evolved GPDs retain a logarithmic ^ de- 
pendence, so that one can also determine a shrinkage 
parameter a' at different scales fi. We find that the 
values of a' for the gluon and quark singlet distribu- 
tions remain close (but not equal) to each other under 
evolution in a model where they coincide at fiQ. In 
an alternative model, where a' for the quark singlet 
is much larger than for gluons at /xq, evolution brings 
their values closer to each other, but clear differences 
remain even at /i^ = 50 GeV^ . An analogous behavior 
is found for /(^) at given ^ in both models. We there- 
fore conclude that one may not take it for granted that 
the t dependence of gluon and sea quark distributions 
is the same at moderate scales. In the flavor non- 
singlet sector, we find that /(^) and a' remain quite 
stable under evolution of the scale. 
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